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Abstract 
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1 INTRODUCTION 



With the production of 10^ Z bosons at LEP rare Z decays with branching ratios as low as 
10~^ come into experimental reach. Among these are decays into heavy quark bound states, 

Z r?^^^Lj + anything, (1) 

where the usual spectroscopic notation is used. Particularly interesting is the possibility to 
produce bottom-charm bound states, so-called B^. mesons. These states have not yet been 
discovered. They are the only quarkonium states with open heavy flavour quantum numbers, 
since the top quark is expected to be too heavy to form and resonances before it decays via 
t bW^. Moreover, Be mesons constitute an excellent laboratory for studying QCD bound 
state dynamics, weak decay mechanisms, quark and gluon fragmentation, and hadronic matrix 
elements. All this provides strong motivation to investigate the prospects for observing Be 
mesons in present and future experiments. 

In the following we briefly describe some characteristic properties of Be mesons, discuss 
the dominant production mechanisms in e~^e~, 77, and pp/pp collisions, and indicate the 
prospects for detection of these new bound states. 

2 SPECTROSCOPY 

The energy levels and other spectroscopic properties of the [be] system can be calculated in 
nonrelativistic potential models or with the help of QCD sum rules. Since the Be resonances lie 
in the mass range between charmonium and bottomonium, the predictions in a given approach 
are constrained by the latter systems. Thus, measurements of Be properties provide interesting 
tests of potential models and sum rule methods. 
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Table 1. Masses and decay constants ||l| of the ^S'o and '^S'l ground states from potential models 
(Martin p[, BT ||^) and moment sum rules [^. The normalization is such that f-^ = 132 MeV 
and fs* = fs^ for nonrelativistic S"— waves. 

In table 1 we exemplify typical predictions for the pseudoscalar (1^5*0) and vector (l^Si) 
ground states. While the mass predictions are very similar in different models and approaches, 
the decay constants are more distinctive (for reviews see e.g. [H). Since the bound state wave 
functions are known in nonrelativistic approximation, one can also calculate more complicated 
quantities such as transition form factors and fragmentation functions, for example, for Be — > 
J/ip and b BeC. We will return to this possibility in section 4. 

3 WEAK DECAYS 

The weak decays of the lightest Be meson proceed via three distinct mechanisms: decay of 
the bottom constituent, decay of the charm constituent, and be annihilation. In fact, the 
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free quark lifetimes, Tb ~ 1.5 and Tc ~ 0.7 ps, are quite comparable, and also the lifetime 
Tftc ~ 5.2 ps due to the annihilation process is not very much longer. The interplay of these 
mechanisms should make decays rather unique and give rise to a rich and interesting decay 
pattern. 

The decay widths corresponding to the above mechanisms have been estimated in quark 
models [^, |^ and using QCD sum rules ^ |[. For the relative rates one finds 

„ „ „_ J 37 : 45 : 18 potentials , . 

Deviations from the ratios 29 : 62 : 8 of the above inverse lifetimes t^^ : t^^ : t^^ are 
expected due to differences in the phase space available for the three decay channels. The 
total Be lifetime is estimated to lie in the range 

= 0.4-0.9ps. (3) 

Decay modes with a. J/ip resonance in the final state provide a particularly favourable 
signature for detection. Unfortunately, the branching ratios are only known with considerable 
uncertainties. For the most important channels one expects 

r 2^5% X = lvi 
Br{Be J/i)X) = \ 0.2 ^ 0.4% X = tt (4) 
{ 0.6^1.1% X = p. 

The inclusive branching ratio for B^. J/ipX is about 20%. It should also be noted that the 
excited states below the open flavour threshold, + m/j ~ 7.142 GeV, will all decay into 
the Be ground state via emission of photons, pions, etc., and hence contribute to the inclusive 
Be rate. 



4 PRODUCTION IN e+e" ANNIHILATION AND 
HEAVY QUARK FRAGMENTATION 

In e~^e~ annihilation B^ mesons emerge from the simultaneous production of a bb and cc pair 
of quarks by strong binding of the be pair in a colour-singlet state. (If no confusion can 
arise, both B^ and B* are simply denoted by B^.) The lowest-order Feynman diagrams of 
this process are shown in Fig. 1. Because of the nonrelativistic nature of \bc] bound states, 
the relative momentum p of the b and c quarks is expected to be small in comparison to the 
quark masses. It is then reasonable to expand the amplitudes in p, and to keep only the 
lowest nonvanishing term |^. In this approximation, the amplitudes for the production of 
5- waves reduce to hard scattering amplitudes multiplied by the 5- wave function ^1^(0) at the 
origin. Similarly, for P-waves one obtains a product involving the derivative of the P-wave 
function. These wave functions can be estimated from potential models. For 5'-waves one can 



also substitute ^'(O) by the corresponding decay constants, |\E'(0)| = ^^JmB^/l2fB^, and use for 
the latter the QCD sum rules estimates indicated in Table 1. 

Straightforward calculation of the Feynman diagrams of Fig. 1 |]10| then yields the total B^ 



and B* cross sections exhibited in Fig. 2. As can be seen, below the Z peak the cross sections 
amount to a few femtobarns only. At the resonance they reach about 1 -r2pb corresponding 
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to 300 ^ 500 events for 10^ Z bosons produced. Above the resonance the cross sections 
again drop rapidly below the level of 1 fb at LEP200. The production process considered 
here dominates other possible mechanisms such as soft fragmentation, that is nonperturbative 
creation of a heavy quark pair from the vacuum, and ly-radiation from light fermions followed 
by W* —>-bc—>- Be- The calculation can be drastically simplified for the following two reasons: 
Firstly, the minimum virtuality k"^ of the gluons in Fig. 1 is determined by the mass of the 
secondary quark pair, i.e. fc^ > 4m^ in la and b, and fc^ > 4m^ in Ic and d. Therefore, the 
gauge invariant subset of diagrams, la and b, is expected to yield the dominant contribution. 
Indeed, at a/s > SOGeV the diagrams Ic and d contribute less than 5% as can be seen from 
Fig. 3. Secondly, in a planar gauge such as Q^G^ = 0, where Q'^ = {pe+ + Pe-Y ^'^d 
denotes the gluon field, the contribution from diagram lb is suppressed relative to that from 
diagram la. Again, this is demonstrated in Fig. 3. Consequently, for s ^ m\ the complete 
process can be described to a very good approximation as production and fragmentation of h 



quarks, e^e —>■ bb, b ^ B^c. This feature has been exploited in [11, 12, 13 




Fig. 1. The lowest-order Feynman diagrams Fig. 2. The total cross sections for e~^e 

, _ _ (*) - 

contributing to e'*~e — > Bcbc. Be be versus the cm. energy. 



Following [|12| , we write the differential B^ cross section as a convolution of the correspond- 



ing bb cross section and a fragmentation function: 

dosXE) = dzda^^i (^^^y") ^^5.(2;,/^), (5) 

where /i is the factorization scale and z = 2EB^I\fs. For [i = 0{mBj the fragmentation 
function can be calculated directly in QCD perturbation theory employing the nonrelativistic 
approximation sketched earlier. However, in order to avoid large corrections proportional to 
In^ in dai,i, it is necessary to take /i = 0{E) and to evolve the fragmentation function Di^^^ 
from jjL ~ niB^ to /i ~ using Altarelli-Parisi equations. To leading order, one has 



din fi'^ 
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with Prh(x) 
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5(1 -x) 



(6) 
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Fig. 3. Cross section ratios calculated from 
Figs. la~d (a), la and h (cfab), O'^d la alone 
in planar gauge (cTa). 



diagrams [o , 
diagram (a' 
D(z) 




Fig. 4. Energy distributions calculated 
from diagrams la-d (solid line), and from di- 
agram la alone in planar gauge (dashed line) 
in comparison to the fragmentation functions 
D{z) given in (8), normalized to unity. 



Fig. 4 shows the fragmentation functions at /i = for b —y Be and b ^ B* normahzed to 
unity in comparison to the distributions as calculated from the complete set of diagrams 
in Fig. 1, and from Fig. la only. With the exception of the threshold behaviour at 2; — > 0, 
the direct calculation of D-^_^^(*) provides an excellent approximation. Evolution to higher 



scales softens these distributions. This is exemplified in |jT^ for ji = lOm^^. However, since 
Jq dz Pii{z) = 0, the integrated fragmentation functions are not affected by evolution. Con- 
sequently, the branching fractions for b B^*^ c are universal and scale independent numbers 



2.2 ■ 10~^ (5c) 
3.1 ■ 10"^ [B*: 



(7) 



where ^^(m^^) = 0.2 has been used. The fragmentation functions for P-waves have been 
worked out in |]T3[. 

It is instructive to compare the perturbative result: 

zil - zY 

D- (,){z,ii) = N M{fi,me)- -F M{z,d), 



v4 



(1 - nz) 



FeAz^d) = —(l + 2d + 2d'^ + 2d^ + d^)-—(2 + 3d+l0d^ + 9d^) 
16 12 

z^ 

+—(1 + 3d+ 17rf2) - 3zd + 1, 
6 
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FBsiz,d) 



z^ 

— (5 + lOd + + 2d^ + d^) (8 + 7ci + d^) 

16 ^4^ ^ 



+ — {9-3d + 3d^) 



z{2 + d) + l, 



(8) 



where d = {nib — ^c)/{^b + ^c), and = mb/{m{, + nic), with the phenomenological 
parametrization introduced in [jl5| for D and i? meson fragmentation functions: 

z(l-z? 



Here, the shape parameter e is considered to be of order m'^/mQ, q and Q denoting the hght 
and heavy constituent quarks, respectively. Interestingly, for e ~ ml/ Ami parametrization 
also provides a quite reasonable description of the distributions In addition to the 
effective shape parameter e, perturbation theory predicts the absolute normalization. 



iV^(*) (/i, mc 



81m? 



(10) 



and the dependence on spin and orbital angular momentum [11, 12, 13, 14 



5 EVENT RATES AT THE Z POLE 

In the fragmentation picture it is easy to estimate the total number of Be mesons produced 
at the Z pole. One simply has to multiply the branching fraction BR{Z bb) = 0.152 with 
the fragmentation probability of b quarks into Be mesons. Since the excited states below the 
BD-threshold eventually decay into Be mesons, they have also to be counted. Adding up the 
S'-waves, one obtains |1^ 

BR{b BecX) ^ 9 ■ 10"^ (11) 

One should, however, bear in mind the considerable uncertainties in this number. These are 
mainly due to the parameters appearing in the normalization factor (0), i. e. the scale /i in 
as, the quark masses and the bound state wave functions or, equivalently, the decay constants. 
Combined, there is certainly a factor two uncertainty in (|TTp as can be seen from the varying 



predictions in the literature ITT], |T2|, |T3|. From the branching fractions given above one expects 
about 

2700 BjBe per 10^ Z. (12) 

About 20% of these decay into final states containing a J/ip. In turn, 12% of the J/tp 
resonances decay leptonically into e'^e" or fi^fi^. Hence, inclusively, one can expect roughly 
60 events per 10"^ Z in the channel Be J/tpX; J/ip ^ l^l^ . It should not be too difficult 
to observe such final states. Unfortunately, the Be signature will be obscured by a large 
background of J /ip events from other sources. Most importantly, one expects more than 3000 
events per 10^ Z from Z — >■ 66; b/b — > J/ijjX; J/ip — > where we have assumed an 

inclusive branching ratio of 1% for b J/ip X and the same for b. This clearly shows that in 
order to discriminate the Be events from the background one has to tag the charm quark jet 
accompanying the \bc] bound states in Z bb; 6 — > ^^^^Ljc; "^^^^Lj — » BeX; Be J/ipX. 

Alternatively, one may search for the Be meson in the particularly clean detection channels 
Be — »• J /ipl'^vi and Be —>■ J /ipTi^ . Using the range of branching fractions given in (|^), one can 
expect 6 -i- 16 events for 10^ Z in the channel and a few events in the channels 

and l^l^p^. These rates are small, but maybe just sufficient for Be discovery. 
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6 Be AT HADRON COLLIDERS 



Given the abundant production rates of 6-quarks in hadronic collisions, the prospects for 
detecting the Be mesons at the TEVATRON or LHC should look much better. In addition to 
the small fragmentation probabilities for b BcC, these rates also afford small Be branching 
ratios. The by far dominant subprocess is gluon-gluon fusion, gg —>■ Bccb, with a small 
contribution coming also from qq annihilation, qq —>■ BcCb. 

Unlike in the case of e~^e~ annihilation, only a minority of diagrams possesses the topology 
of fragmentation processes, e. g. gg ^ bb]b —>■ B^c. Correspondingly, heavy quark fragmenta- 
tion can only be expected to dominate Be production at large transverse momenta. At smaller 
Pt and for total cross sections, recombination processes such as gg bb cc; be Be will play 
an important role. The onset of the asymptotic regime in p^, where the fragmentation picture 
is valid, is a priori not clear. We will return to this point in the next section. 







(Tdj,.(cuts) 


A^i„c(cuts) 


TEVATRON 
{25 pb-^) 


5.3 nb 


PT > 10 Gey, |y| < 1 
0.13nb 1.6-10^ 


LHC 
(100/6-1) 


60 nb 


PT > 20GeV, \y\ < 2.5 
0.5 nb 2.1-10^ 



Table 2. Cross sections and number of events at ^/s = 1.8 TeV (TEVATRON) and 14(16) TeV 
(LHC) for inclusive (inc) and direct (dir) Be production. 

In Table 2 we quote estimates for cross sections and production rates which can be found in 
the literature. The total cross sections aine have been calculated in [|TB] taking into account the 
complete sets of 0(a^)-Feynman diagrams for the gg and qq subprocesses, and summing the 
individual cross sections for IS* and 25* waves. On the other hand, the high-p^- cross sections 
Cdiricuts) have been obtained in [|1^] using an approximation analogous to (|^), and considering 
only direct fragmentation into the Be ground state described by (||). Also applied are usual 
rapidity cuts. We emphasize that the uncertainties in these predictions are even larger than 
in the case of e^e" annihilation. The main reason is the scale ambiguity in the QCD coupling 
constant which enters here in fourth power. This makes the evaluation of higher-order QCD 
corrections very desirable. 

As anticipated, even at high pt the production rates are large enough to afford detection 
channels with branching ratios of order 1% and below. In fact, CDF has already recorded some 
candidate events [jl8|. A B-physics experiment at LHC should certainly be able to investigate 



the Be system in some detail. 



7 PHOTON-PHOTON PRODUCTION - A THEORET- 
ICAL STUDY CASE 

As pointed out, hadronic production of Be mesons involves fragmentation and recombination 
mechanisms and is thus more involved than Be production by e^e~ annihilation. In order to 
learn about the relative importance of fragmentation and recombination and their distinctive 
features, we have investigated Be production in 77 — >■ Bebc. This process is simpler than gg 
fusion, but involves the same two mechanisms as can be seen from the Feynman diagrams in 
Fig. 5. Moreover, it is interesting to find out whether there is a chance to observe Be mesons 



in laser induced 77 collisions feasible at future linear collider facilities [19 
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In Fig. 6 we compare the cross section derived from the complete setQ of 0{a'^al) diagrams 
indicated in Fig. 5 with the result obtained from the gauge invariant subset (I), which includes 
the fragmentation diagrams. Also shown is the cross section resulting from a convolution of 
(t(77 — > bb) with the fragmentation function -Dj^Bc analogy to (^). We observe that the 
integrated cross section is dominated by the recombination mechanism. This is not unex- 
pected. We also find agreement between the predictions of diagrams (I) and simple 6— quark 
fragmentation. Apparently, there is a gauge in which the first diagram of set (I) dominates. 




Fig. 5. Different topologies of the lowest- Fig. 6. The total cross section for 77 — > Bcbc 

order Feynman diagrams contributing to 77 as function of the cm. energy. 

Bed. 

Furthermore, we have found that the energy distribution resulting from subset (I) is very 
similar to the distribution predicted by the pure fragmentation description. This is illustrated 
in Fig. 7 together with the complete energy distribution including also the recombination 
processes (II). The latter completely change the character of the energy spectrum giving rise 
to rather soft Be production. 

We have also compared the corresponding prp distributions. As expected, the recombination 
processes mainly populate the low-p^ region. Nevertheless, their pj- spectrum has a long tail 
towards high pt- For example, at ^/s = 100 GeV and px = 20 GeV, the contribution from 
the diagrams (I) and (II) are still equally important. Recombination dies out completely only 



at Pt ~ AO GeV. These findings put some doubt on the estimates in |17| of Be production at 
high-pT in hadronic collisions, where the fragmentation approximation is used at px <^ \/s. 

In contrast to the energy and px distributions, the angular distributions resulting from 
fragmentation and recombination processes are quite similar in shape. All distributions are 
sharply peaked for | cos6'| — > 1, with the distribution derived from the complete set of diagrams 
in Fig. 5 being somewhat fiatter than the approximations from diagrams (I) and from da^i ® 

As a final remark, for ^/s > 100 GeV we predict an integrated cross section below 12 fb. 



■^Diagrams of type (I) with the photons coupled to the c-quark line are not taken into account. They 
contribute an extra 15% to the total cross section iBOl 
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This shows that 77 coUiders cannot be expected to provide promising opportunities in Be 
physics. 



da 
adz 
3.0 



2.0 



1.0 



0.0 

0.0 0.2 0.4 0.6 0.8 2 1-° 

Fig. 7. Energy distributions calculated from the diagrams (I) and (II) of Fig. 5, and from the subset 
(I) alone in comparison to the fragmentation function D^_^q^{z, fj,) given in (8), normalized to unity. 



8 SUMMARY 

The Be system provides plenty of possibihties to investigate interesting aspects of strong and 
weak interactions. Most importantly, the quarkonium nature of these bound states allows 
to calculate and study also genuinely nonperturbative quantities such as total cross sections, 
decay amplitudes, transition form factors and decay constants. Therefore, Be physics can 
serve as a testing ground for quantitative theoretical approaches to confinement problems. 
Practically, Be mesons are coming into experimental reach, may be at LEP, likely at the 
TEVATRON, but most certainly at LHC. 
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